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ABSTRACT 

We present new Chandra X-ray and Giant Meterwave Radio Telescope (GMRT) radio observations of the 
nearby merging galaxy cluster Abell 754. Our X-ray data confirm the presence of a shock front by obtaining 
the first direct measurement of a gas temperature jump across the X-ray brightness edge previously seen in the 
imaging data. A 754 is only the fourth galaxy cluster with confirmed merger shock fronts, and it has the weakest 
shock of those, with a Mach number M= 1 .57^ [|. In our new GMRT observation at 330 MHz, we find that the 
previously-known centrally located radio halo extends eastward to the position of the shock. The X-ray shock 
front also coincides with the position of a radio relic previously observed at 74 MHz. The radio spectrum of the 
post-shock region, using our radio data and the earlier results at 74 MHz and 1 .4 GHz, is very steep. We argue 
that acceleration of electrons at the shock front directly from thermal to ultrarelativistic energies is problematic 
due to energy arguments, while reacceleration of preexisting relativistic electrons is more plausible. 
Subject headings: galaxies: clusters: general — galaxies: clusters: individual (A754) — intergalactic medium 
— radio continuum: galaxies — X-rays: galaxies: clusters 



1. INTRODUCTION 

Galaxy clusters form via mergers of smaller substructures. 
During such mergers, most of the kinetic energy of the gas 
belonging to the colliding subclusters is dissipated by shocks 
and turbulence into thermal energy of the intracluster medium 
(ICM) of the resulting system. Shocks and turbulence are 
also expected to amplify the cluster magnetic fields and ac- 
celerate cosmic ray particles from the thermal ICM, or reac- 
celerate preexisting relativistic particles. These non-thermal 
components manifest themselves as diffuse synchrotron radio 
sources, known as radio halos and relics (e.g., Ferrari et al. 
2008; Cassano 2009 for recent reviews), and inverse Comp- 
ton X-ray emission at high energies (e.g., Fusco-Femiano et 
al. 2004; Rephaeli & Gruber 2002; see, however, Wik et al. 
2009). 

Shock fronts represent a unique observational tool to study 
the physical processes in the ICM. They create sharp discon- 
tinuities in the cluster X-ray surface brightness images and al- 
low to measure the gas velocities in the sky plane using X-ray 
imaging spectroscopy (e.g., Markevitch, Sarazin, & Vikhlinin 
1999). Shock-heated regions are routinely observed in merg- 
ing clusters. However, known shock fronts are still rare, be- 
cause one has to catch the front when it has not yet moved 
to the outer, low surface brightness regions of the cluster, 
where the background X-ray emission dominates. Moreover, 
the merger has to be occurring nearly in the plane of the sky, 
otherwise projection effects could hide the gas density and 
temperture jumps. 

Until now, reliable detections of merger shock fronts have 
been reported in only three galaxy clusters. One is the Bul- 
let cluster (1E0657-56, Markevitch et al. 2002), another is 
Abell 520 (Markevitch et al. 2005), and the two fronts have 
been recently discovered in Abell 2146 (Russell et al. 2010). 
Only in these clusters both the sharp gas density edges and 
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the unambiguous temperature jumps were found, allowing the 
identification of the brightness feature as a shock and deter- 
mination of the gas velocity. In this paper, we report on the 
new shock front in the merging cluster Abell 754. We present 
the analysis of Chandra observations and study the connec- 
tion between the shock and the diffuse cluster radio emission, 
through Giant Metrewave Radio Telescope {GMRT) observa- 
tions at 330 MHz and Very Large Array (VIA) archival data 
at 1.4 GHz. 

We assume a flat cosmology with Ho = 70 km s -1 Mpc" 1 and 
n Q = 0.3, in which 1" is 1 .054 kpc at the redshift of A754. We 
adopt the convention S oc v~ a for the synchrotron spectrum. 
Uncertainties are 68%, unless stated otherwise. 

2. THE MERGING CLUSTER A754 

A754 is a rich nearby cluster at z = 0.0542 (Struble & Rood 
1999) in the stage of a violent merger. It has been actively 
studied in the optical and X-ray bands, and is considered the 
prototype of a major cluster merger. Previous studies revealed 
that the cluster has a complex galaxy distribution (Fabricant 
et al. 1986, Zabludoff & Zaritzki 1995), X-ray morphology, 
and gas temperature structure (Henry & Briel 1995; Hen- 
rikesen & Markevitch 1996; Markevitch et al. 2003, here- 
after M03; Henry et al. 2004). These data indicate that A754 
is undergoing a major merger of two components along an 
east-west axis, probably with a non-zero impact parameter. 
Hydrodynamic simulations were able to reproduce most of 
the observed X-ray features in A754 by considering an off- 
axis collision between two subclusters with a mass ratio of 
2.5:1 (Roettiger et al. 1998). Chandra data suggested that the 
merger may be more complex, possibly involving a third sub- 
cluster or a cloud of cool gas decoupled from its former host 
subcluster (M03). 

Krivonos et al. (2003) reported an edge-like surface bright- 
ness feature in the ROSAT PSPC image of the cluster, located 
east of the core, which looked like a shock front. They mod- 
eled it with a radial density profile with a jump, deriving a 
Mach number M = 1 .71^ 24 fr° m this jump under the assump- 
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tion that this is indeed a shock front. This brightness feature 
was also apparent in the Chandra image by M03, but it was 
at the edge of the Chandra field of view, which did not allow 
its detailed study. Henry et al. (2004) presented a tempera- 
ture map of the cluster from an XMM-Newton observation. 
The presumed post-shock region was shown to be hot, as ex- 
pected for a shock front. However, because the X-ray bright- 
ness ahead of this front is very low, the XMM observation did 
not have sufficient sensitivity to determine the temperature in 
the pre-shock region, in order to detect the expected tempera- 
ture jump. They were only able to derive a temperature in the 
annulus around the cluster that included the pre-shock region 
(region 12 in Fig. 6 of Henry et al. 2004), but that annulus 
is dominated by emission from regions on the other side of 
the cluster, unrelated to the shock front. We note that many 
features that look like shocks in the X-ray images turned out 
to be "cold fronts" (Markevitch & Vikhlinin 2007), and the 
difference between these phenomena is the sign of the tem- 
perature jump at the feature. In this work, we present a direct 
temperature measurement in the pre-shock region from a new 
Chandra observation. 

A754 also exhibits evidence for ultrarelativistic electrons 
and magnetic fields that coexist with the ICM. Fusco-Femiano 
et al. (2003) reported a 3er hard X-ray excess at E > 45 
keV with BeppoSAX, which may come from inverse Compton 
scattering of the Cosmic Microwave Background (CMB) pho- 
tons on relativistic electrons. A clear evidence of nonthermal 
emission in A754 comes from the radio observations. A radio 
halo was detected at 74 MHz, 330 MHz and 1.4 GHz with the 
VIA (Kassim et al. 2001, hereafter K01; Bacchi et al. 2003, 
hereafter B03, and references therein). Moreover, the pres- 
ence of two peripheral radio relics were reported in K01, East 
and West of the radio halo. Only the eastern emission was 
seen in B03. The whole cluster diffuse radio emission has 
been recently studied in the frequency range 150-1360 MHz 
by Kale et al. (2009). 

3. CHANDRA OBSERVATIONS AND THE SHOCK 
FRONT 

A754 was first observed with Chandra ACIS-I in October 
1999 (OBSID 577). This pointing was centered on the cluster 
center and used by M03 to study the cluster temperature struc- 
ture. In Febraury 2009, a new long Chandra ACIS-I observa- 
tion (OBSID 10743) was pointed on the putative shock front, 
in order to derive the temperature profile across the shock. 
The new exposure partially overlaps the old one (see Fig. |T|l, 
and in this paper we use both observations for the spectral 
analysis of regions where they overlap. 

Both datasets were processed and cleaned in a standard 
manner (e.g., Vikhlinin et al. 2005). The data were cleaned 
of flares and the blank-sky datasets were normalized as de- 
scribed in Markevitch et al. (2003b). OBSID 577 had small 
flare contamination and the final clean exposure was 39 ks. 
No background flares were present in the new observation, 
so we used the full exposure of 95 ks. The old observation 
was taken in FAINT mode, while the new one in VFAINT, 
which allowed additional background filtering to be applied. 
The detector + sky background was modeled using the blank- 
sky background datasets corresponding to the dates of ob- 
servations (periods B and E, respectively). We normalized 
the blank-sky backgrounds using the the ratio of counts in 
the high-energy band 9.5-12 keV, which is free from the sky 
emission. This correction was within 10% of the exposure 
ratios, as expected. 



All point sources were masked out for the extraction of 
brightness profiles and spectra. The instrument responses for 
spectral analysis were generated weighting the detector ARF 
and RMF with the cluster brightness within each spectral ex- 
traction region (Vikhlinin et al. 2005). We used the most re- 
cent calibration products — namely, version N0008 for the 
telescope effective area, N0006 for the CCD quantum effi- 
ciency and gain files, and N0005 for the ACIS time-dependent 
low-energy contaminant model. We also assessed the system- 
atic uncertainty of our results by using a newer, experimental 
update to the ACIS contaminant model (A. Vikhlinin, private 
communication; see below). 

Figure[T]shows a slightly smoothed ACIS image of the clus- 
ter in the 0.5^4 keV band, obtained by co-adding the two 
observations and correcting for exposure nonuniformity. As 
found by previous studies, the cluster undergoes a complex 
merger with the main axis along the NW-SE direction. The 
image clearly shows a brightness edge to the east of the dense 
elongated core, perpendicular to the merger direction. This 
is the putative shock front reported by Krivonos et al. (2003). 
The new observation provides sufficient statistics to determine 
the exact nature of this feature. 

3.1. Density profile across the edge 

Figure [2] shows a radial X-ray surface brightness profile 
across this edge, extracted in a 25° sector shown in Fig.[T] To 
facilitate proper geometric modeling, the sector is centered on 
the center of curvature of the brightness edge (dashed white 
circle in Fig. [T]l and encompasses its most prominent segment 
(avoiding the apparent decrease of the edge brightness con- 
trast at greater opening angles). All point sources are care- 
fully masked. The energy band for the profile was restricted 
to 0.5^1 keV to minimize the dependence of X-ray emissivity 
on temperature and maximize the signal-to-noise ratio. Only 
OBSID 10743 is used for the brightness profile, since the ear- 
lier observation does not cover the pre-shock region. The ex- 
posure of this observation alone is more than sufficient for 
the post-shock region. We consider only the region outside 
the bright elongated core, which is obviously unrelated to the 
shock. 

The brightness profile across the edge (Fig. [2^) exhibits the 
typical shape of a projected spherical density discontinuity 
(Markevitch et al. 2000), and we fit it with such a model, un- 
der the assumption that the curvature along the line of sight 
is the same as in the image plane. The model radial density 
profile that we use consists of two power laws, p oc r a ' and 
p oc r° 2 , on two sides of the edge, with an abrupt jump at the 
edge (Fig.|2|3). For each set of the profile parameters, we pro- 
jected the corresponding emission measure profile onto the 
image plane and fit to the observed X-ray brightness profile 
(Fig. |2^), treating the slopes and the radius and amplitude of 
the jump as free parameters. The fit was restricted to the in- 
teresting radial range — staying as close as possible to the 
edge on the inside in order to avoid contamination from un- 
related core structure, but extending as far as possible on the 
outside for proper deprojection (see dashed ticks on the outer 
sides of the sector in Fig. [TJ. As seen from Fig. [2^, the model 
fits very well. The best-fit density jump, after a 3% correc- 
tion of the 0.5 - 4 keV Chandra brightness fo r the measured 
temperature difference across the edge (see § 3.2 below), is 



Pi/ Pi = l-80^Jj f 5 . The confidence interval is evaluated by al- 
lowing all other model parameters to be free. 

3.2. Temperature profile across the edge 
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FIG. 1. — Chandra slightly smoothed 0.5^4 keV image of A754 from a sum of the new 95 ks expo sure and t he archival 39 ks observation. Also shown is the 
25° sector used to extract the X-ray surface brightness profile and for the spectral analysis (see § |3.1| and § |3.2| ). Dashed ticks outside the sector mark the radial 
range used for the profile fitting; small ticks inside the sector mark the radial bins chosen to derive the temperature profile; bigger red ticks indicate the best-fit 
density jump, measured from its center of curvature (marked by the dashed circle). 



TABLE 1 
X-RAY TEMPERATURE FITS 



Region obsid 577 obsid 10743 Simult. fit 

Current (N0005) ACIS contaminant model: 

1 17.fi*-* 9.0«° 10.13-g 

2 9.(,-;. 14.2+|f 12.5+J! 

3 8.7^-2 16.4+11 13.3lf| 

5 4.5^ 

3deproj .... ... ... l6A -3.l 

Experimental ACIS contaminant model: 
3 8.7+2| 14.2+2-8 12.2+2-' 

4 ' 



We now extract spectra from several radia l bins in the same 
sector across the density edge (Fig{T[ §3.1) to determine the 



gas temperature jump across the edge. There is enough sta- 
tistical accuracy to divide the low-brightness region outside 
the edge into two bins, and the inner brighter region into three 
bins, avoiding getting to close to the cluster cool elongated 
core. The five radial bins are marked by ticks inside the sec- 
tor (FigJTJi; the longer red ticks mark the best-fit position of 
the density jump (Fig{T[ The first high-brightness re- 

gion adjacent to the jump extends across the best-fit jump po- 
sition a bit, in order to avoid contaminating the adjacent low- 
brightness region by any irregularities in the edge shape. The 
spectra were fit in the 0.8-9 keV band with XSPEC, using the 
absorbed thermal plasma model WABS(APEC), with elemen- 
tal abundances as free parameters, fixing the absorption col- 
umn to the Galactic value (N H = 4.82 x 10 20 cm" 2 , Kalberla et 
al. 2005). We also tried to free Nh, and all fits were consistent 
with the Galactic value. All fits had acceptable \ 2 values. The 
systematic uncertainty of background modeling was evaluated 
by varying the blank-sky background normalization by ±2% 
(68%, Hickox & M arkevitch 2005) and added in quadrature to 
the statistical uncertainties. It was negligible for the brighter 
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FIG. 2. — (a) The 0.5-4 keV surface brightness profile across the front in the sector shown in Fig.^(see text). Errors are lcr; histogram shows the best-fit 
model that corresponds to the radial density profile shown in panel (b). The dashed line marks the best-fit position of the shock front (measured from its center 
of curvature). 
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FIG. 3. — Temperature profile across the shock (errors are 68%). Solid 
crosses show projected temperature fits; dashed line for the first post-shock 
region shows a deprojected temperature (modeling the projected contribution 
of the pre-shock gas based on the best-fit density model). For post-shock 
regions, simultaneous fits for two observations are shown, while two outer 
regions are covered only by the recent observation (see Table [T] for more 
details). Dotted vertical line shows the position of the shock front from the 
density fit (see Fig.[2j. 



regions inside the edge, but similar to the statistical uncer- 
tainty for the two low-brightness bins (the observation expo- 
sure was selected to achieve this). 

Solid crosses in Fig. [3] show the resulting temperature val- 
ues, using the simultaneous fit to both observations for the 
inner regions and OBSID 10743 for the outer regions. The 
figure shows a clear temperature jump of the sign that corre- 
sponds to a shock front. 

Table[T]gives separate fits for two pointings; the regions are 



numbered in order of increasing radius in Fig. [3] We note 
that two Chandra observations are in mild disagreement for 
the post-shock regions (though none deviates by more than 
2cr from the simultaneous fit). This may be caused by resid- 
ual calibration problems — these regions are at the edge of 
the ACIS-I field of view in OBSID 577, while at the center in 
OBSID 10743. In addition, these observations are separated 
by 10 years. We tried an experimental update to the recently 
released time dependence of the time-dependent ACIS-I con- 
taminant model (A. Vikhlinin, private communication) to see 
if it changes the results qualitatively, for the two bins adjacent 
to the front. The results are shown in Table [T] The disagree- 
ment is slightly reduced. Most importantly, the temperature 
jump does not change qualitatively, and is also seen in the 
recent observation alone, using either calibration. To the ex- 
tent the results from different spatial regions can be compared, 
our post-shock high temperatures are in broad agreement with 
those derived with XMM by Henry et al. (2004). 

Thus, the temperature profile confirms that the brightness 
edge is a shock front. From the Rankine-Hugoniot jump con- 
ditions, we can use our best-fit gas density jump to derive a 
Mach number of this shock. Assuming monoatomic gas with 
7 = 5/3, we obtain M = 1 .57^;}| 

The temperature jump can give an independent Mach num- 
ber estimate, though usually with a lower accuracy. To check 
for the consistency of these estimates, we first get a depro- 
jected temperature in the first post-shock radial bin, since in 
the spherical geometry, the temperature we measure there is 
affected by projection of the cooler emission from the outer, 
pre-shock bins. (At the same time, the pre-shock tempera- 
ture profile appears consistent with isothermal, so the true 
temperature should be close to the projected one, to a suffi- 
cient accuracy). From the best-fit density model, we estimate 
a projected emission measure fraction in bin 3 from the outer 
regions 4 and 5 in Fig.[3]to be 21% and 3%, respectively. Ne- 
glecting the latter and fitting the spectrum for region 3 using 
an additional thermal model with the normalization and tem- 
perature (including its uncertainties) corresponding to region 
4, we obtain a "deprojected" temperature shown by the dashed 
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cross in Fig. [3] 

The ratio of the post-shock to pre-shock temperatures is 
thus 2.2^o', which corresponds to M = 2.1^ \ (68%), consis- 
tent (within its expected larger uncertainties) with the value 
from the density jump obtained above. 

4. RADIO OBSERVATIONS 

In this section, we present our GMRT observations of A754 
at 330 MHz, and analysis of the archival VLA data at 1.4 GHz. 
Details of observations are summarized in Table [2] which re- 
ports the telescope, project code, observing date, frequency, 
total bandwidth, total time on source, synthesized half-power 
bandwidth (HPBW) of the full array, rms level (Ict) at full res- 
olution, u—v range, and the largest detectable structure (LDS). 

4.1. GMRT observations at 325 MHz 

A754 was observed with GMRT at 325 MHz in June 2005 
for a total time on source of 2.5 hr (Table EJ. The observa- 
tions were carried out using the upper and lower side bands 
simultaneously (USB and LSB, respectively), for a total ob- 
serving bandwidth of 32 MHz. The default spectral-line ob- 
serving mode was used, with 128 channels for each band and 
a spectral resolution of 125 kHz/channel. The LSB dataset 
was corrupted and could not be processed, so only the USB 
data, which are centered at approximately 330 MHz, were 
used to produce the images presented here. The dataset was 
calibrated and analyzed using the NRAO Astronomical Image 
Processing System package (AIPS). We refer to Giacintucci 
et al. (2008) for a complete description of the data reduction 
procedure used in this paper. 

Due to the large field of view of GMRT at 330 MHz (pri- 
mary beam ~ 1 .8°), we used the wide-field imaging technique 
at each step of the phase self-calibration process, to account 
for the non-planar nature of the sky. We covered a field of 
~ 2.7° x 2.7° with 25 facets. The final images were produced 
using the multi-scale CLEAN implemented in the AIPS task 
IMAGR, which results in better imaging of extended sources 
compared to the traditional CLEAN (e.g., Clarke & Ensslin 
2006; for a detailed discussion, see Appendix A in Greisen, 
Spekkens, & van Moorsel 2009). We used three circular 
Gaussians as model components. One of the Gaussian was 
chosen to have zero width to accurately model point sources 
and small-scale structures, and the other two have a width of 
20" and 45" respectively, to progressively highlight the ex- 
tended emission during the clean. 

Beyond the image at full resolution (10.0" x 9.1"), we pro- 
duced images with lower resolution (down to ~ 100"), taper- 
ing the u — v data by means of the parameters robust and 
uvtaper in the task IMAGR. Even though only half of the 
data was usable, the sensitivity of the final images is quite 
good — the rms noise level (la) ranges from 0.45 mJy beam" 1 
in the full resolution image to ~1 mJy beam" 1 in the lowest 
resolution images. We estimate that the flux density calibra- 
tion uncertainties are within 5%. 

4.2. VLA archive data at 1.4 GHz 

VLA observations of A754 at 1 .4 GHz (project AF372) were 
presented by B03. We extracted these observations from the 
archive and re-analyzed them in order to ensure the best possi- 
ble comparison with our GMRT data. The observations were 
obtained using the D configuration and two IFs, centered at 
1365 MHz and 1435 MHz (see Table |2]for details on the ob- 
servations). Standard calibration and imaging were carried 



out using AIPS. The dataset was self-calibrated in phase only. 
The higher frequency IF was found to be affected by strong 
radio interference, which compromises the quality of the im- 
ages produced using both IFs. For this reason, only the low- 
frequency IF was used for the analysis presented in this paper. 
The final images were obtained implementing the multi-scale 



clean option in IMAGR, as for the 330 MHz data (see f 4.1 



The rms noise in the image at full resolution is 50 /iJy beam 
(Tabled. A slightly higher noise (la ~ 60-70 /iJy beam" 1 
was achieved in the low-resolution (HPBW~70") images. 
The average residual amplitude errors in the data are of the 
order of < 5%. 

5. RADIO ANALYSIS 

In Fig. [?] we present the new^\GMRT image at 330 MHz 
at the resolution of 63" x 52" (left panel) and the VLA 1365 
MHz image at full resolution (HPBW=64" x 39"), overlaid 
on the optical POSS-2 red frame of A 754 (right panel). The 
first positive contour corresponds to the 3a level in both im- 
ages (i.e., 3.3 mJy beam" 1 at 330 MHz and 0.21 mJy beam" 1 
at 1365 MHz). Numerous discrete radio sources are located in 
the cluster area. Some of these sources are embedded within 
the diffuse emission. We identify all the radio galaxies found 
by B03. Following their same notation, these are labelled 
from SI to S8 in both images. The point source S5 is un- 
detected at 330 MHz. 

In order to properly image the cluster diffuse emission, we 
subtracted all the discrete radio galaxies in the cluster region 
from the u — v data at both frequencies. We first subtracted 
the brightest radio sources in the field (including the extended 
radio galaxy north-east of SI in Fig. [4]), to simplify the detec- 
tion and subtraction of the fainter discrete galaxies embedded 
in the diffuse emission (from SI to S8). 
We finally produced a set of low-resolution images from the 
"subtracted" data sets at both frequencies, to highlight the ex- 
tended emission only. 

Our results are shown in Fig. [5] The left panel shows 
contours of our VLA low -resolution image at 1365 MHz 
(HPBW=70" x 70"), before the subtraction of discrete radio 
sources^] The rms noise of the image is 60 /iJy beam" 1 (1 
a); this is significantly lower with respect to the noise level 
reached by B03 in their image of same resolution (see Fig. 
3 of B03). In the right panel of Fig. [5] we present images of 
the diffuse emission only, obtained after the subtraction of dis- 
crete sources described above. The GMRT 330 MHz intensity 
contours at the resolution of 109" x 74" are overlaid on the 
VLA 1365 MHz image (grayscale), produced with a restoring 
beam of 70" x 70". Both images show the very large extent 
and complexity of the diffuse emission in A754. The diffuse 
source spans ~ 1 .4 Mpc in the east-west direction and ~ 750 
kpc along the north-south axis, as measured from the 3a con- 
tour level at 330 MHz. The 1365 MHz image reveal a ridge in 
the radio brightness, extended westbound the 330 MHz con- 
tours. On the other hand, the high frequency emission is less 
extended in the eastern region (see below). We point out that 
these are reliable features, not due to a shift between the two 
images. Indeed we carefully checked the astrometry of our 

3 These observations, available from the GMRT public archive, has been 
also recently used by Kale et al. 2009, for their low frequency study of A 754. 
However, they do not show the radio images obtained from these data. 

4 Here S9 ? indicate the peak of emission visible to the north-east of S7 
in Fig. [4] This apparently coincides with an optical counterpart, and might be 
an additional point source; however we considered it as a peak of the diffuse 
emission and we did not subtract it from the data. 
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TABLE 2 

Summary of the radio observations. 





Radio 
telescope 


Project 
code 


Observation 
date 


V 

(MHz) 


Au 
(MHz) 


t 

(min) 


HPBW, p.a. 

("x", °) 


rms 
(/Jy beam -1 ) 


u-v range 
(kA) 


LDS 

(') 


GMRT 
VLA-D 


08GBA01 
AF372 


June 23, 2005 
Sept. 25, 2000 


325 
1365/1435** 


32(16)* 
50 


150 
160 


10.0x9.1,-64 
64.5x38.8,5 


450 
50 


~ 0.08-25 
~ 0.13-4.7 


~32 
~ 15 



Notes to Table|2] * Only the USB dataset was used for the analysis (see Sect. |4.l| . ** Only the the lower frequency IF was used for the analysis (see Sect. |4.2| 
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FIG. 4.— Left: GMRT 330 MHz image of A754 at the resolution of 63.5" X 52.1", p.a. 30°. The lcr level is 1.1 mJy beam" 1 . Contours start at -3.3 mjy 
beam -1 (dashed) and +3.3 mjy beam -1 , and then scale by a factor of 2. Right: VLA full-resolution image at 1365 MHz (contours), overlaid on the optical POSS-2 
red image of A754. The restoring beam is 64.5" X 38.8", p.a. 5°, and the lcr noise level is 50 /xjy beam -1 . Contours are spaced by a factor of 2, starting from 
±0.15 mjy beam -1 (negative contours are shown as dashed). In both images, labels indicate the discrete radio galaxies identified at 1.4 GHz by B03. 
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FIG. 5. — Left: VLA low-resolution image at 1365 MHz. The restoring beam is 70" X 70", p.a. 0; contours are spaced by a factor 2 starting from ±3<r, the 
lcr level is 60 fjjy beam -1 . Labels indicate the discrete radio galaxies identified within the diffuse emission. Negative contours are shown as dashed. Ticks mark 
contours that correspond to holes in the radio brightness. Right: diffuse emission after the subtraction of all identifiable discrete sources. Red contours show a 
low-resolution GMRT image at 330 MHz with a restoring beam of 109" X 74", p.a. 43° and the lcr noise level of 1 mjy beam -1 . Contours are spaced by a 
factor of 2, starting from 3<r. Greyscale shows a low-resolution VLA image at 1365 MHz — same as in left panel but without the discrete sources. The peak flux 
density is 3.77 mjy beam -1 and the lowest level of greyscale corresponds to ~ 2<r. The green ellipse is the area covered by the 2cr contour (0.4 Jy beam -1 ) of the 
radio relic in the VLA image at 74 MHz (K01). Ticks mark contours that correspond to holes in the radio brightness. 
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FIG. 6. — GMRT low resolution radio contours (same as in Fig. [3] right panel), overlaid on the Chandra image of A754 (same as Fig.^ see §[5J. The dashed 
contours show the 3<r negative values; the red arc indicates the bestlit position of the density jump (see § |3.1) . 
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radio images at both frequencies, by comparing the positions 
of the peaks of the radio sources with that of the optical coun- 
terparts. 

The green ellipse in Fig. [5] (right panel) represents the 2a 
contour (0.4 Jy beam" 1 ) of the eastern radio relic found in the 
74 MHz VIA image by KOI. The 330 MHz diffuse emission 
extends eastward to the outer side of this relic . Interestingly, 
we notice that the shock front found in § |3.2| coincides with 
the edge of the emission at 330 MHz (see Fig. |6l), possibly 
suggesting a connection between the radio edge and the shock 
(see §[6) . 

We detect significant emission at 1365 MHz only in part of 
the region covered by the green ellipse (see Pig . J5]> . Our VIA 
images at 1365 MHz show more diffuse emission, south of 
the S2-S4 sources, with respect to B03. This is most likely 
due to differences in the editing of th e data, in the imaging 
procedure (multi-scale clean, see § |4.2| i and in the accuracy of 
the subtraction of discrete sources (see also the comment in 
B03); our final rms value is indeed a factor 1.6 better. 
We do not detect any radio emission in the region of the west- 
ern relic found by KOI. 

5.1. Spectrum of the diffuse radio emission 

The radio images in Fig. [5] (right panel) show a complex, 
Mpc-scale diffuse source associated with A754. We use our 
330 MHz GMRT and the reprocessed 1365 MHz VIA data 
presented in 341 along with the VIA image at 74 MHz kindly 
provided by NT Kassim (see Fig. 1 in KOI), to derive its in- 
tegrated spectrum. The flux densities of the diffuse emis- 
sion, after the subtraction of the discrete sources, are 828±41 
mJy and 89±5 mJy at 330 MHz and 1365 MHz, respectively. 
These are also reported in Table [3] along with the associ- 
ated uncertainties and angular resolution of the images used 
for the measurements. We obtained these values by integrat- 
ing over the same area that encompasses the whole diffuse 
emission as detected at both frequencies [^] We note that the 
flux density at 1365 MHz is lower than the value reported 
in B03 where, however, the authors notice that the subtrac- 
tion of discrete sources is not accurate (B03; see also § |5j. 
The flux density of the whole diffuse emission at 74 MHz is 
~ 6.6 Jy (Tab. 13). This value was integrated over the same 
area described above and does not include the contribution of 
the brightest radio galaxies at the cluster center (i.e., S2, S3 
and S4, and the extended radio galaxy north-east of SI; see 
Fig. H) which amounts to a total of ~ 2.9 Jy (as estimated in 
K01 ^ However, due to the lower sensitivity, only part of the 
extended emission visible at 330 MHz and 1365 MHz is de- 
tected at 74 MHz. Thus, our estimate should be considered as 
a lower limit to the flux. The spectral index for the whole dif- 
fuse emission between 330 and 1365 MHz is a = 1.57 ±0.05. 
This value is steeper than typically found for giant radio halos 
(a ~ 1.2—1.3; e.g., Ferrari et al. 2008), although a number 
of halos with ultra-steep spectrum (a > 1.5) has been discov- 
ered in the past few years (e.g., Bmnetti et al. 2008, Brentjens 
2008, Macario et al. 2010). 

For simplicity, in the following we will call relic the feature 
discovered by K01 at 74 MHz, that they defined as eastern ra- 
dio relic (see also discussion, §[6). 

In our GMRT image at 330 MHz, the relic does not appear 
as a distinct feature separated from the radio halo, nor a clear 
increment of the surface brightness of the diffuse emission is 

5 We also checked that no significant change in the spectral index occurs 
by integrating within the emission at the 2<r level. 



detected in that region. This may imply that the spatial seg- 
regation between the relic and the halo in the 74 MHz image 
(see Fig. 1 in K01) may be due to a steeper spectrum of the 
relic at v < 300 MHz combined with the lower sensitivity of 
the 74 MHz observation. 

To obtain the radio spectrum of the relic region, we integrated 
the flux densities at 330 MHz and 1365 MHz over the area 
defined by the 2a contour level of the relic at 74 MHz (green 
ellipse in Fig. [5] right panel). We obtained 106±5 mJy at 330 
MHz and 6.0±0.3 mJy at 1365 MHz (see Table |3J. 

In Fig. [7] we show the integrated radio spectrum of the 
relic, in the frequency range 74-1365 MHz. The relic has in- 
deed a very steep spectrum, with oi^q^, = 1.77 ±0.10, and 
a o = 2-02 ± 0.04. An extrapolation of this power law to 
lower frequencies predicts a 74 MHz flux density that is con- 
sistent with that reported in K01. For a direct comparison, the 
flux densities of the whole diffuse emission are also shown. 

A spectral study of the diffuse emission in A 754 was pre- 
sented in Kale et al. (2009). However their analysis is limited 
to specific regions of emission detected in their GMRT 150 
MHz image (that do not show the whole extended emission). 
Therefore, we cannot compare our spectral results to those 
found by them. 



TABLE 3 

Flux densities of the diffuse radio emission in A 754 



v (MHz) HPBW , 



S„ (mJy) 



Ref. 



relic 



whole emission 



74 
330 
1365 



316 X 232 
109 X 74 
70.0 x 70.0 



1489±223 
106±5 
6.0±0.3 



~6600 
828±41 
89±5 



this work; K01 
this work; Fig. 5 
this work; Fig. 5 



10" 



whole diffuse emission 



10 r 



1UU T r.„, n 

Log v [MHz] 

FIG. 7. — Integrated spectrum of relic region (rilled dots) in A 754 in the 
frequency range 74-1365 MHz. The solid line connects the values of flux 
density reported in Table[3] The spectrum of the whole diffuse emission be- 
tween 330 and 1365 MHz is also shown (open triangles and dashed line); the 
estimated flux density at 74 MHz (filled triangle) is in reasonable agreement 
with the extrapolation of the spectrum. 
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6. DISCUSSION 

The main result of the present paper is the clear X-ray de- 
tection of a shock front in the intracluster gas of A 754, which 
coincides with the edge of non-thermal radio diffuse emission 
at 325 MHz. 

Particle acceleration at the shock front cannot be respon- 
sible for the whole Mpc-scale diffuse emissi on s een in Fig. 
[5] Given the shock parameters derived in § |3.2| the down- 
stream velocity of the gas is ~ 1100 km s -1 , which implies 
that the size of the diffuse emission produced by the shock 
can be only ~ 100-200 kpc, once the radiative lifetime of the 
emitting electrons (1 - 2 x 10 8 yr) is taken into account (see 
also the discussion in Markevitch et al. 2005 and Brunetti et 
al. 2008). However, the spatial coincidence of the ultra-steep 
spectrum radio edge with the shock front suggests a physical 
connection between these two features, and possibly an indi- 
rect connection between the shock and the diffuse emission 
on the larger scale. If one assumes that relativistic electrons 
in the relic region are accelerated via Fermi mechanism b y th e 
shock front with M ~ 1.6 (as observed in X-rays, see §3.2| i, 
the expected slope of the radio spectrum from that region, ac- 
counting for the downstream losses, is a ~ 2.3 (a = a,> !7 + 1 /2, 
where a inJ = (p-l)/2 and p = 2(M 2 + 1)/(M 2 - 1); e.g., Bland- 
ford& Eichler 1987 ), in rough agreement with the radio ob- 
servation (see §5.1| l. However, we can rule out direct shock 
acceleration on energy grounds. If we extrapolate the power 
law spectrum of electrons responsible for the relic emission 
down to thermal energies (as we should if the seed electrons 
are thermal), the energy in relativistic electrons will be of the 
same order as the thermal energy in the ICM requiring an im- 
plausible high acceleration efficiency. This energy argument 
implies that observable radio emission cannot be produced at 
shocks with M < 2 (Hoeft & Brriiggen 2007). 

Below we discuss possible scenarios for the origin of the 
cluster-scale diffuse emission in A754. 

1 . The edge emission might arise from relativistic plasma 
that is re-energized at the shock. If fossil electrons are 
present in the ICM, compression by the shock can sig- 
nificantly increase their synchrotron emission at the ob- 
serving frequency (e.g., Ensslin et al. 1998, Ensslin & 
Gopal-Rrishna 2001, Markevitch et al 2005). Adiabatic 
compression increases the maximum synchrotron fre- 
quency emitted by a population of fossil electrons by 
r 4 / 3 times (e.g., Markevitch et al. 2005), where r is the 
shock compression factor (considering a shock moving 
perpendicular to the line of sight). In our case, this im- 
plies only a moderate boosting of a factor 2. The result- 
ing synchrotron emission may light up at our observing 
frequency. The drawback of this scenario is that in this 
case the fossil plasma must have been injected in the 
ICM only a few 10 s yrs ago, because the fossil elec- 
trons must exist at energies just below those necessary 
to emit at the observed frequency, and at those energies 
the lifetime is short. 

This problem remains even if one assumes that the fos- 
sil plasma is confined (not mixed) by the ICM (Ensslin 
& Gopal-Krishna 2001), due to the low Mach number 
of the shock. 

2. Shock reacceleration of fossil relativistic plasma can 
potentially be an efficient process even in the case of 
weak shocks, though this process at weak shocks is 
poorly understood and its efficiency is uncertain. As- 



suming the standard linear shock acceleration theory 
(e.g., Blandford & Eichler 1987), the spectrum of reac- 
celerated electrons is: 

N(j) = (p+2)j- p f N fss (x)x"- l dx (1) 

where p = 2(M 2 + 1)/(M 2 - 1) is the slope of the spec- 
trum of accelerated electrons immediately after the 
shock, Nf S s is the spectrum of fossil electrons in the 
upstream region, and 7 ml >, is the minimum energy of 
the fossil electrons at which the reacceleration process 
works. If the fossil relativistic electrons have an expo- 
nential cutoff in their power spectrum (as they should if 
they result from radiative cooling, e.g. Sarazin 2002), 
shock reacceleration should create a spectrum above 
that cutoff with a slope similar to the one in the clas- 
sical Fermi reacceleration. Ass uming M ~ 1.6, consis- 
tent with our findings in § |3.1| a power-law tail of ra- 
dio emission with a ~ 2.3 can be generated at energies 
7 > 1000 (Figure [§}. This scenario avoids the previ- 
ously mentioned energy argument, because the result- 
ing steep spectrum of electrons does not extend down 
to the thermal energies. 




Log(7) 



FIG. 8. — Spectrum of reaccelerated electrons according to Eq.[T]We also 
correct for the shock compression factor in the downstream region, r. The 
thick solid line is the initial spectrum of fossil electrons (assuming an age 
of the population of a few Gyrs and typical cluster physical parameters, e.g. 
Sarazin 2002). Different models show the spectrum for Mach numbers 1.3 
(dotted line), 1.57 (dashed line), 2.5 (long-dashed line), and 3.5 (dot-dashed 
line). A minimum energy 7,,,,,, = 10 is adopted in the calculations. 

3. The connection between the edge and the shock might 
be indirect. For instance, the passage of the shock may 
have driven small-scale turbulence in the ICM, which 
may be long-lived and also produce the radio halo on 
the larger scale. In this case, the radio edge would 
mark the region where turbulent acceleration is just be- 
ginning to occur. This should happen after ~1 eddy 
turnover-time of the turbulence, w L/V/ (where L and 
Vi are the maximum scale and velocity of the turbulent 
eddies), since the turbulent modes at smaller scales are 
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the most important in the particle acceleration process 
(e.g., Brunetti & Lazarian 2007 and references therein). 
This may also explain the very steep spectrum of the 
radio emission in the shock region. A possible indirect 
connection between shocks and the cluster-wide diffuse 
radio emission is suggested by the present radio obser- 
vation, as well as those of several other clusters, where 
bridges of emission connecting radio relics and halos in 
several clusters are seen (see, for instance, the discus- 
sion in Brunetti et al. 2008). A theoretical exploration 
of this scenario deserves a future effort. 

7. CONCLUSIONS 

We have presented a combined X-ray and radio study of 
the nearby merging galaxy cluster Abell 754. 
The new Chandra observation confirms the presence of a 
merger shock front east of the cluster core by providing 
the Mach number from the density jump across the shock 
M= 1.57±g;[2, and a direct measurement of a gas tempera- 
ture jump of the right sign across the arc-like brightness edge, 

r 2 /r 1 =2.2!';i. 

The new GMRT radio image of the cluster at 330 MHz re- 
veals that the centrally located diffuse emission is very ex- 
tended and complex, and is mainly elongated in the east- 
west direction. Most interestingly, the eastern edge of this 
source coincides with the shock. We studied the spectral prop- 
erties of the diffuse radio emission, using the GMRT data, 
archival VIA data at 1.4 GHz, and previous VIA results at 
74 MHz. We find that the region next to the radio edge (the 
relic region in K01) has a very steep integrated spectrum, 

with ao3Gft = 2m ± °- 04 and a omcHz = ini ± °- 10 - T^s is 
steeper than the average spectrum of the whole diffuse emis- 
sion, cxq ^chI = 1-57 ± 0.05, which is in itself steeper than the 
typical spectrum of radio halos. 



The spatial coincidence of the steep spectrum radio edge 
and the merger shock front suggests a physical connection. 
The low Mach number of the shock and the very steep syn- 
chrotron spectrum of the edge allow us to rule out the scenario 
of direct shock acceleration of thermal particles for the origin 
of the edge, because it would require an implausible high ac- 
celeration efficiency (of order 1). 

Possible alternative scenarios are reacceleration and/or adia- 
batic compression of fossil relativistic plasma in the ICM. Of 
these, shock reacceleration is more plausible for this weak 
shock. It should produce a spectrum of electrons that is in 
agreement with the observed steep radio spectrum of the relic 
region. 

Finally, an indirect connection between the shock and the ra- 
dio edge is also possible. The shock passage may drive small- 
scale turbulence that may reaccelerate electrons over the clus- 
ter volume. 

Upcoming GMRT deep low-frequency observations will al- 
low us to perform a more detailed study of the diffuse emis- 
sion and help to determine the most plausible scenario for its 
origin. 
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